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TURBULENT -SKIN- FRICTION AND HEAT-TRANSFER CHARTS 


ADAPTED FROM THE SPALDING AND CHI METHOD 

By Luther Neal, Jr., and Mitchel H. Bertram 
Langley Research Center 

SUMMARY 

Charts are presented which allow a rapid determination of local and average skin 

friction and heat transfer on flat plates in air. The charts cover a Mach number range 
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from 0 to 20, a Reynolds number range from 10 to 10 in decade increments of the expo- 
nent, and a wall -temperature — stagnation-temperature ratio range of 0.1 to the adiabatic- 
wall case. 


INTRODUCTION 

In the design of large hypersonic vehicles, the need for accurate estimates of turbu- 
lent skin friction and heat transfer is important. A number of prediction methods have 
accumulated over the years, but the results from these methods often differ widely and 
recourse to experiments must be made. In early 1964, Spalding and Chi (ref. 1) compared 
experimental skin friction with 20 different turbulent theories. A theory of Van Driest 
(ref. 2), employing the Karman mixing length, was found to yield the lowest root-mean- 
square error (11.0 percent). However, Spalding and Chi were able to formulate a method 
applicable to a flat plate that gave an error of 9.9 percent. In a paper by Bertram and 
Neal (ref. 3), a body of heat- transfer data in the Mach number range from about 5 to 9 and 
wall-temperature — stagnation-temperature ratios down to about 0.1 was found to support 
the predicted effect of wall -temperature — stagnation-temperature ratio from the Spalding - 
Chi method. Because of current interest in the Spalding-Chi method, it was decided to 
present, in convenient graphical form, the results of extensive calculations that were 
accumulated while preparing reference 3. 

SYMBOLS 

C F average skin -friction coefficient based on conditions at edge of boundary layer 

Cf local skin-friction coefficient based on conditions at edge of boundary layer 



F c> F Rx’ F R 0 functions in the Spalding-Chi theory given by equations (3), (2), and (1), 
respectively 

M Mach number at edge of boundary layer 

Np r Prandtl number 

Ngt local Stanton number based on conditions at edge of boundary layer 

Rjj Reynolds number based on distance x and conditions at edge of boundary 

layer 

T absolute temperature 

t = T/T e 

Subscripts: 

aw adiabatic wall 

i incompressible value 

w wall 

t total 

e edge of boundary layer 


DISCUSSION 

In developing their method, Spalding and Chi postulated that a unique relation exists 
between F c Cj and FpR where Cf is the local skin-friction coefficient, R is the 
Reynolds number, and F c and Fp depend only on Mach number and wall-temperature — 
stream -temperature ratio. The function F c was calculated by means of the mixing- 
length theory, and the function F R was deduced semi-empirically. Thus, values of skin 
friction could be obtained but heat transfer was not considered. 

Tables and charts for F c and Fp to be used in the calculation of skin friction 
were presented; however, for extensive use of the theory, it was found more convenient to 
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compute values of F c and Fr. When the Reynolds number is based on momentum 
thickness, F R = Frq and is given by the empirical expression (eq. (52) of ref. 1). 


. -1.474. 0.772 

* R0 - % l aw 


( 1 ) 


or when based on length of turbulent flow x, Fj^ = Fr x and is given by equation (25) of 
reference 1. 


Frx = 


f R0 

Fc 


( 2 ) 


Only the Reynolds number based on distance x is considered in this paper. To compute 
F c , equation (39) of reference 1 was integrated to yield the expression 


F c 



2^-aw " l(*aw *w “ M 2 ~ ^aw “ *wj] 
(taw + tw)^ ” 4t w 
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( 3 ) 


The proper quadrant in which to evaluate the inverse sine is shown in the following sketch: 



M 


In the limit as M approaches 0, equation (3) reduces to 



By using these expressions and the variation of F c Cf and F c Cf with FrxR x , 
as given in table 7 of reference 1, and extended where necessary by equations (28) to (30) 
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of reference 1*, values of Cf and C F were computed. Stanton numbers were calcu- 
lated by the modified Spalding-Chi method described in appendix A of reference 3. This 
approach consists of applying the Karman form of the Reynolds analogy in the incom- 
pressible plane so that the Karman equation is modified to be 


2N S t 

Cf 



(5) 


When the Karman factor was applied in this paper, Np r was taken as 0.725. 

Although directly applicable to a flat plate, the results for local skin friction and 
heat transfer may be applied to a cone with a known or assumed location of transition by 
using the transformation given in appendix B of reference 3. Results given in reference 4 
(pp. 25 to 30) may aid. in the application of the results to obtain the average skin friction 
on a cone or tapered wing. 

PRESENTATION OF RESULTS 

Charts of Cf/Cf^, Ngt/Ng^j, and Cjr/Cp^ are presented in figure 1 as a func- 
tion of Mach number for R x = 10^ to 10^ and Tw/Tt = 0.1 to the adiabatic- wall value. 
The recovery factor was taken to be 0.89. The charts are normalized by using the incom- 
pressible values of skin- friction coefficient given by the Karman- Schoenherr formulas: 


= l°gl 0 (c F)i Rx) = log 10 (2Ri 

C 0.242C F ,i 

1,1 0.242 + 0.8686jc^T7 


^Equation (29) of reference 1 actually contains two typographical errors; the correct 
form which was determined by referring back to the paper in which the expression was 
originally presented is 

R x,i = ^("g) 4 + (k 3 e) _1 {[6 - 4(kuq) +(kug)^6xp(Ku^) - 6 - 2(ku*) 

- IlK ) 4 - U Ku $ 5 - e^K ) 6 - M Ku $} 

+ /2 V l /2 

where Uq = (— j , K = 0.4, and E = 12. 
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The normalizing values of Ng^ ^ are those obtained from the unmodified Karman equa- 
tion (eq. (A4) of ref. 3) in conjunction with equation (6) and are presented in figure 2. 

In figure 1 for the adiabatic- wall cases at M = 0, the value of unity is not neces- 
sarily obtained because the Karman-Schoenherr formulas were used for normalizing 
rather than the incompressible values of Spalding and Chi. Actually, the incompressible 
values of skin-friction coefficients of Spalding and Chi do not differ significantly from 
those of Karman-Schoenherr as seen in figures 3 and 4, but because of the widespread 
use of the Karman-Schoenherr formulas in normalizing experimental data in past publica- 
tions, they were chosen for use here. For some cases, the skin-friction or heat -transfer 
ratio drops below the value 0.1 before Mach number 20 is attained. Extrapolation is 
aided by noting that the logarithm of these ratios is essentially linear with the logarithm 
of M between Mach numbers of 10 and 20 at the high Reynolds numbers. Thus, extrap- 
olation may be simply done on log-log paper. 

It should be noted that at the higher Mach numbers and lowest Reynolds numbers, 
the equivalent incompressible Reynolds numbers are much lower than those expected for 
turbulent flow to exist. (See ref. 5.) 

CONCLUDING REMARKS 

Utilizing the Spalding-Chi turbulent -skin-friction theory, charts are presented which 

allow a rapid determination of the local and average skin friction on flat plates in air. 

An extension was made to allow the determination of heat transfer by modifying Spalding- 

Chi’ s skin -friction results according to the Karman form of the Reynolds analogy factor. 
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The charts cover a Mach number range from 0 to 20, a Reynolds number range from 10 
9 

to 10 in decade increments of the exponent, and a wall-temperature — stagnation- 
temperature ratio range of 0.1 to the adiabatic -wall case. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 13, 1967, 

129-01-09-14-23. 
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Figure 2 - Variation of incompressible Stanton number with Reynolds number based on the Karman form of the Reynolds analogy equation 

and the Karman-Schoenherr skin-friction equation. Np r = 0.725. 
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